iV-dimensional nonlinear Fokker-Planck equation with time-dependent coefficients 
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An TV- dimensional nonlinear Fokker-Planck equation is investigated here by considering the time 
dependence of the coefficients, where drift-controlled and source terms are present. We exhibit the 
exact solution based on the generalized gaussian function related to the Tsallis statistics. Further- 
more, we show that a rich class of diffusive processes, including normal and anomalous ones, can be 
obtained by changing the time dependence of the coefficients. 
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Anomalous diffusion processes appear in a large class of 
systems in several contexts. Some illustrative examples 
are diffusion in plasmas jlj, relative diffusion in turbulent 
mediaQ, CTAB micelles dissolved in salted water 0, sur- 
face growth and transport of fluid in porous media|4| , two 
dimensional rotating flowQ, subrecoil laser cooling^, 
diffusion on fractals 0, anomalous diffusion at liquid 
surfaces Q, diffusion in linear shear flows fl , enhanced 
diffusion in active intracellular transport fH|7particle dif- 
fusion in a quasi-two-dimensional bacterial bathjll , and 



dence of the diffusion coefficient [J17LJ1 
poral dependence of the drift termll 




and tem- 



spatiotemporal scaling of solar surface flows [Q, among 
others. 

The existence of the anomalous diffusion and its ubiq- 
uity has motivated, in particular, the analytical study 
by consideringnonlinear[jl3[ ^ [H], [jj, [0| and 

fractionalEO, Fokker-Planck equations, spatial depen- 



In this brief report, we consider an iV-dimensional non- 
linear Fokker-Planck equation incorporating the time de- 
pendence in every coefficient, where drift-controlled (ex- 
ternal force) and source terms are present. In our exact 
solution, a rich class of anomalous behaviors can arise by 
choosing appropriate time dependence of the coefficients. 
These results indicate that the possible anomalies in dif- 
fusive processes can appear as a consequence of different 
causes. This paper contains, essentially, the main results 
obtained in Refs. |ll| |l8) and |24| as special cases. 

In order to investigate, in an exact way, a large class of 
anomalous diffusions, let us consider the A^-dimensional 
nonlinear Fokker-Planck equation 



J^(r, t) = V(t)V 2 [p(r, t)Y - V • [F(r, t)p(r, t)} - a(t)p(r, t), 



(1) 
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where we incorporate the time dependence in the external 
force, 



F(r,i) =ki(t) + fe 2 (t)r, 



(2) 



in the diffusion coefficient T>(t), and in the source term 
a(t). In particular, we consider the case where the space 
constant term in F(r,i) can be taken with different co- 
efficients. This situation is useful, for example, to study 
diffusion in the gravitational field [ki = (0, 0,ki z )]. The 
harmonic potential is considered isotropic. 

The source term in Eq.(|l|) can be removed by an ap- 
propriate change in the solution: 



p(r, t) = cxp 



- / a(t')dt' 



(3) 



This way, p(r, t) obeys the equation 

j t p(r, t) = D (t)V 2 [p(r, t)Y V • [F(r, t)p(r, t)] , (4) 

with D(t) = V{t) exp [(1 - v) f* a(t')dt'] . Thus, Eq. (g) 

has the same structure of Eq. ([!]) without the source 
term, but with an additional time dependence of the dif- 
fusion coefficient. Observe that this additional time de- 
pendence of D(t) is induced by the nonlinear term p", 
disappearing when v = 1. 

In order to obtain an exact solution for Eq. (Q) with 
the external force term given by Eq. (||), we are going to 
employ the ansatz 



P(r,t) = 



Z(t) 



1 - (1 - <7)/3(i) (r-ro(i)r 



V(i-?) 



(5) 



2 



if 1 - (1 - q)P(t)(r - r ) 2 > 0, and p(r,t) 



if 



It is usual to identify normal diffusion process by a lin- 



1 — (1 — q)P(t)(r — r ) 2 < (cut-off condition). We ear growth in time of the variance a 2 =< (r — ro) 2 >. 



would like to remark that Eq. (|5|) can be justified via 
dimensional analysis and related to the renormalization 
group theory [^5). Furthermore, this ansatz can also be 
connected with Tsallis statistics^, £27|. In addition, Eq. 
(||) reduces to the Gaussian when q — ► 1. In fact, by 
defining the function exp q (— x 2 ) = [1 — (1 — q)x 2 ] 1 ^ 1 ^ qS> 
if l-(l-q)x 2 > 0, and exp g (-a; 2 ) = if \-{l-q)x 2 < 
as a q-gaussian, we obtain the usual gaussian function by 
taking the limit q — ► 1. Eq. (^) is a solution of Eq. (Q) 
when v — 2 — q and the time dependence of (3{t), Z(t), 
and ro(i) = J2iLi x 0i(t) e i are ruled by the following sys- 
tem of equations: 



1 dZ 
Z~dt 
1 dp 

P~dt 



2(2 - q)NDf3Z- 1+q - Nk 2 , 
-4(2-q)DZ- 1+q (3 + 2k 2 , 



and 



dxoi 
dt 



= ku - k 2 x 0l 



(6) 
(7) 

(8) 



Note that Eqs. (||) and (^) are nonlinear, Eq. (||) is 
independent of /3(i) and Z(t), and ku(t) does not appear 
in the nonlinear coupled differential equations for (3{t) 
and Z(i). Furthermore, the spatial independent term in 
the external force only affects the time dependence of xqi. 
So, Eq. (§) leads to 



x i(t) 



2^(0) + / k u (s)e-^ds 



(9) 



where p(t) = f Q k 2 (s)ds. For example, in a 3-dimensional 
space with the presence of an isotropic time-independent 
harmonic potential and the gravitational field, i.e., 
k 2 (t) = k 2 = constant, k\ x — k\ y — and k\ z {t) = k\ = 
constant, we get x {t) — x o (0)e~ k2t , yo{t) = y (0)e~ k2t , 
and 



zo(t) 



MO) + | (e 



k 2 t 



1 



-kit 



(10) 



The solution for the nonlinear coupled equations for 
Z(t) and (3(t), Eqs. (|) and (|), is given by 



and 



Z(/. = z 



N/ci 



-2/ci 



(11) 



(12) 



with d = 2 + N(l - q), /3 = (3(t = 0), Z = Z(t = 0) 
and 



f{t) = e - Cl ^ l) x 



iVfc 2 (s)-2iV(2-g)/3 Z 9 - 1 D( S ) 



(13) 



Other time dependences on a 2 are commonly related as 
anomalous diffusion, for instance, supcrdiffusive, subd- 
iffusive, exponentially diffusive, and localized. In our 
study, we obtain a large class of diffusive processes which 
include these examples. In addition, tr 2 can have different 
behaviors for small and large times, enabling the descrip- 
tion of a rich structure of diffusion regime. Here we are 
going to analyze, as an illustration, a set of representa- 
tive kinds of these anomalous behaviors in the asymptotic 
regime taking some specific dependence of the coefficients 
into account. In this direction, from Eqs. (|3j> and (||), 
we investigate the asymptotic temporal behavior of the 
variance 



J(r-r Q ) 2 p(v,t)d N r 
fp(r, t)d N r 



(14) 



where C(q, N) is a constant depending only on q and N. 
Note that the convergence of the integral for q > 1 in a 2 
imposes a restriction over the parameters: 2 + iV(l — q) > 
2{q — 1). This implies that c\ is a positive constant for 
all q values. 

Firstly, let us consider the case without the source 
term, with the diffusion coefficient constant, and with the 
time dependence of the harmonic external force given by 
k 2 (t) = k t~ b . From Eq. Q, we obtain 



/(*) 



-dkt 1 -" /(l-b) 



c 2 e 



-cikt 1 -" i '(l-b) 



c 3 t b 



(15) 



for large t and b < 1, where c 2 and C3 are constants 
which depend on N, fc, q, and b. Again, for large t and 
b = 1, we have /(*) - t- 2k - c 2 t- 2k (t 1+2k - 1)/(1 + 2k), 
and finally, for b > 1, we get /(<) ~ t. By using the 
fact that the mean square displacement is given by a 2 ~ 
/(i) 2 / Cl , several asymptotic behaviors can be obtained. 
We summarize in Tab. (|) the possible behaviors related 
to the above asymptotic results. When we restrict our 
analysis to the one-dimensional linear case (v = 1), this 
drift-controlled anomalous diffusion contains the results 
given in p^] as a particular case. 

Consider now the nonlinear diffusion equation with 
neither source nor linear external force, but with the 
time dependence of the diffusion coefficient given by 
D(t) = V(t) = D t d . In this case, for large t, Eq. (|l|) 
leads to 



f(t)~t d+1 . 



(16) 



Since a 2 ~ /(t) 2 /I 2 + Ar ( 1 -'?)] , there is a competition be- 
tween the parameters q and d to define the diffusion 
regimes. Tab. (||) contains a summary of these regimes. 

Another possible situation is to take the nonlinear 
diffusion equation with a time-dependent source term, 
a(t) = aot a , with time-independent diffusion constant, 
T>(t) = Dq, and without linear external force, k 2 (t) = 0. 
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TABLE I: Large time behavior of a 2 ~ f{t) 2/ci for a(t) = 0, 
D(t) = V(t) = D , and k 2 {t) = fct -6 , where ci = 2 + AT(1 - 
q) > 0. 
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t 2/cl 


Ci-diffusive a 


6 = 


k>0 


(1 _ g-cifct-j2/ C1 


Ornstein-Uhlenbeck 


6 = 


k < 


e l\k\t 


exponentially diffusive 


< 6 < 1 


fc > 


t 2b/c 1 


ci-diffusive 6 


< & < 1 


A; < 


e 2\k\t'-» 


less than 








exponentially diffusive 


6 < 


fc > 




localized 


6 < 


fc < 




more than 








exponentially diffusive 


6= 1 


k > -1/2 


t 2/cl 


ci-diffusive a 


6=1 


fc = -1/2 


(ilni) 2/ci 


log divergent 


6= 1 


k < -1/2 


f 2|fc| 


superdiffusive 


b> 1 


k 




ci -diffusive" 



"The process is superdiffusive for ci < 2, normal for ci = 2, and 
subdiffusive for ci > 2. 

6 The process is superdiffusive for ci < 2b, normal for ci = 2b, and 
subdiffusive for ci > 26. 



TABLE II: Large time behavior of a 2 ~ f(t) 2/ci for a(t) = 0, 



0. 



V(t) = Do* , and k 2 {t) = 0, where ci = 2+N(l-q) > 



d 
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a 2 (t) 


description 


d = 


<1 


= l 


t 


normal 


d = 


'I 


< l 


t 2/l2+JV(l- 9 )J 


subdiffusive 


d = 


'I 


> l 




superdiffusive 


d < 


<1 


= l 




subdiffusive 


d > 


q 


= l 


t i+d 


superdiffusive 


d>N(l-q)/2 




q 


t s(i+rf)/ia+JV(i-9)j 


superdiffusive 


d<N(l-q)/2 




q 


t 2{l+d)/[2+N(l-q)\ 


subdiffusive 



For v = 1 (q = 1) the source term does not affect the 
diffusion regime. On the other hand, for !/ ^ 1 (9 ^ 1) 
and large t, Eq. (|l3|) gives 



/(t)~t°exp 



(g - l)a ^+ Q 
1 + a 



(17) 



when a > -1. For a = -1, we have fit) - t(?-i)«o+i ; 
and for a < -1, Eq. (pj) reduces to /(*) ~ Tab. 
(|H) gives us a summary of the above behaviors. To 
conclude our observations about the anomalous diffusion 
induced by the time-dependent coefficients, we stress that 
the investigation of a more complex time dependence of 
the coefficients can be reduced to the analysis of Eq. ( 13 ) . 

Summing up, we have investigated an ./V-dimensional 
nonlinear Fokker-Planck equation by incorporating the 



time dependence every coefficient, including those of the 
external force and the source term. An exact solution 
is obtained in the case of external force with isotropic 
spatial linear term and a possible anisotropic spatial con- 

TABLE III: Large time behavior of a 2 ~ f(t) 2/ci for a(t) = 
a t a , V(t) = D , and k 2 (t) = 0, where a = 2 + JV(l-g) > 0. 



a 2 (t) 



description 



(q - 1)oq < a = 



,(q-l)a t) t_ 1 
(q-l)a 



— 



stationary 



(q - Vjao > a = 



,(9-l)a *_ 
(q-l)a 



exponentially 
diffusive 



a = -1 



t 2l(g-l) ao -l)/ C1 Cl _diffusive a 



less than 

2[(<;-l) Q0 ]t 1 + °] 

(q — l)ao > — 1 < a < e c i< 1 + a ' exponentially 

diffusive 



(q - 1)qo < -1 < a < 



pTST7ci 



ci -diffusive 



(g - l)ao > a > 



2[(g-l)a ]t 1 + a ] 
g ci(l + a) 



more than 

exponentially 
diffusive 



(q - l)a < a > 



£-2a/ci 



localized 



a < 



— 



ci-diffusive c 



"The process is superdiffusive for c\ < 2[(q — l)ao — 1], normal 
for ci = 2[(g — l)ao — l]i and subdiffusive for c\ > 2[(q — l)ao — 1]. 

'The process is superdiffusive for c\ < 2|o|, normal for ci = 2|a|, 
and subdiffusive for cj > 2|a[. 

c The process is superdiffusive for c\ < 2, normal for c\ = 2, and 
subdiffusive for c\ > 2. 



stant term. This anisotropic constant term is useful, for 
instance, to analyze diffusion when the gravitational field 
is relevant. We have showed that a rich class of anoma- 
lous behaviors arises by choosing appropriate time de- 
pendence of the coefficients. These results indicate that 
the anomaly in diffusive processes may appear as a conse- 
quence of different causes. In particular, the combination 
of nonlinearity and time dependence of the coefficients 
can lead to the normal diffusion (time linear growth of 
the mean square displacement). This fact implies that 
the normal diffusion, in general, can not be associated 
with the Gaussian shape of the distribution p(r,t), i. e., 
the linear time increase of variance does not necessarily 
mean that we are in the presence of ordinary diffusion. 
Finally, we hope that the results obtained here may be 
useful to clarify a possible origin of a large class of dif- 
ferent anomalous diffusive processes in theoretical and 
experimental contexts. 
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